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ABSTRACT 


5 report was prepared to describe the results of 
Enc cperimental study of Їгее turbulent mixing in axial 
En 126 gradients. Both the pressure field and the velocity 
field were extracted using pitot-static probes and both a 
single normal wire and an X-array hot-wire anemometers. 
The resulting hot-wire data shown in Figures 1-80 present 
excellent visual representation of the physical aspects 
er lurbulent flow. 

Саал == оте Ее results of this experiment show 
. πονεῖν that the Reynolds Stresses in the flow are 
ШО ед бу some function to the turbulent kinetic energy; 
however, no specific function was found that would produce 


Eucocrelationship observed in this study. 
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NOMENCLATURE 


reussimesqsapparent Viscosity 
X-Array hot-wire lead (utv) 
X-Array hot-wire lead (u-v) 
Constant in hot-wire calculations 
Constantin apparent viscosity relation 
шосе rbUlent Shear coefficient 
Correlation coefficient 

Wake generator diameter 

Linear output from hot-wire 

Shape factor 

ВЕ От j direction (1=1,2,... 
M pulente kinetic energy 

Prandtl mixing length 

Pressure 

Dynamic pressure 

General flow parameter 

Apparent viscosity constant 
Characteristic time 

Time 

Spatial Coordinates 

х - Axial 

y - Lateral 

Corresponding velocities 


Uhn-cedgctendjstie velocity 


j=1,2,.. 





B Constant пп Hot-wire calculations 


==. 


б Displacement thickness 
6, Energy thickness 
€ Mixing length parameter = КП 
6 Momentum thickness 
и Шаппааг viscosity 
Hr Apparent viscosity 
ν E mamarzskıpetie Viscosity 
Ve Apparent kinetic viscosity 
p Density 
Т Reynold's Stress 
V Differential operator DEL 
E iImteeran from a to b with respect to y 
a 
= Zee zgerıvarıvezwich respect со г 
E Tamary derivative with respect to r 
| | Indicates absolute value 
Pußscripts 
A Parameter of lead A on hot-wire 
B Parameter of lead B on hot-wire 
e Local external flow conditions 
х,у ПОПЕО РЕ a relatar tterentiation with respect to x,y 
Т Shear related quantity 
superscripts 


— 


Indicates a time averaged quantity 


' Designates a fluctuating quantity 





INTRODUCTION 


MOSTE real fluid flows the transition from laminar to 
Eusutent flow carries the fluid mechanist from a situation 
where he can safely predict the values of the flow-field 
parameters from the solution of the Navier-Stokes Equations 
memone where predictions are obtainable only оп the average, 
οὗ Yet most practical applications require at least 
meme consideration of turbulence and how it can be dealt 
With. In the past the engineer has assumed various plausible 
hypotheses in order to introduce a measure of order and 
111 to turbulent flow. It is the intention of 
Mis study, then, to attempt to experimentally verify certain 
of these hypotheses. Before describing the currently 
accepted assumptions regarding turbulent flow, a general 
ἡ οἱ their history will be presented to provide the 
mecessary background information. 

A turbulent flow is defined as one in which an irregular 
ul cbpuablon or mixing motion is superimposed on the princi- 
Ф tluid motion. Its onset is accompanied by a noticeable 
сисе in the law of resistance. Іп the case of laminar 
pipe flow, the longitudinal pressure gradient which main- 
pm (ШЕ погіоп 15 proportional to the first power of the 
Eb Contrast, in turbulent flow this pressure 
gradient becomes nearly proportional to the square of the 


mean flow velocity. In order to calculate flow field values 





Ma Ur robvulent field, then, some logical hypothesis must 
Badge tO account for this change in the apparent viscosity 
of the flow. First, it is generally assumed that the fluid 
motion can be separated into a mean and a Mera tine part 
ШЕП аге linearly related to the total fluid motion in 


such a way that a flow parameter (r) can be defined 
r=r+p T) 


iu e definitions can be substituted into the equations of 
motion to prođuce a new set of equations involving mean and 
Mec cuatineg terms and their cross products. Assuming that 
ши-нопте атегаре of the fluctuations is zero, the equations 
may be averaged to get the following momentum equation for 

a two-dimensional mean flow (see Appendix A for the complete 


development): 


e за du Dex 
—dU — Qu _ e 0 u P 
ua + y еи + Ve (x) ο. (А-22) 


where the apparent kinematic viscosity Су, (х) Io a ως σας 
Миг curbulent flow field, and is calculated by using some 
logical empirical hypothesis το ο ο Ἱπαξτιως Its value. 
Since only the mean values of the flow parameters are 
known, some means of relating these values to the Reynolds 
stresses produced by the fluctuating part of the flow is 


required. The first hypothesis for relating these quantities 





ШӘ ргоровеа by J. Boussinesqg, 


“ 


ze 
ei 


T = ~pu'v' = A 


O τ 


where A. 15 ап apparent or eddy viscosity that is not a 
mea property but depends on the structure of the flow 


Mela. 


A = p v, (x) 


This apparent viscosity is generally assumed to be propor- 


tional to some characteristic velocity (v") times the mean 


mMelocity gradient: 


* 2u Gre) 


Ши 1925 Prandtl made an important advance in turbulent 


theory by introducing the concept of mixing length and the 


gelation, 


(I-4) 


where 
£ = mixing length 


9 IS = characteristic velocity 





Mis relation has produced relatively accurate values for 
, Some ecc ases. however, this definition of the kine- 
matic viscosity implies that it should vanish where the 
pa Pl c Ly gradient is Zero, Шира 15 сопсгааїссогу to experi- 
mental evidence. Several subsequent hypotheses were proposed 
taylor, Von Karman, and Prandtl himself, but Sen were 
инт ша сес 1 application to free turbulent flow. 

ντ ο ορ со develop a universally valid 
πες] Гог the function Hg > ргашес етос тосшоей an energy 
posed hypothesis which was subsequently improved on by 
Emmons (1954), and Glushko (1965), 


1/2 


αρ κ г (1-2) 


Hr 


(ж КО > 


N 
T 


Ec was found to be a better characteristic velocity than 
immune Carlier mixing length hypothesis, but the length 
Scales were difficult to ascertain. More modern theories, 
referred to as two-equation viscosity models, use sophisti- 
K analytical methods to determine this length scale. 

The apparent viscosity is first assumed to be related to the 


kinetic energy. Then some hypothesis of the form 


10 





ШЕ пасе со make it possible to calculate the characteristic 
Busch. One of the most accepted of these relations was 


introduced by Chou in 1945: 


е = 3272 ΠΡ, 
и, = C1 р 222 
Thus 
2 1/2 
шиг 11: 9 (1-7) 


It was this type of assumption which initiated the pro- 
ject reported here. Using both pitot-static and hot-wire 
ENemometer sensors, an attempt was made to validate the 
wo Nesis that shear was predominantly a function of the 


turbulent kinetic energy. 


TL 





11. EQUIPMENT 


A. WIND TUNNEL 
M cceuerad Deseription 

lnescunmnertocomersts essentially of three major 
Euponents: the inlet, the test section, and the blower as 
shown in Figures 81-83. The inlet was designed to contract 
Slow moving inlet air to a high-velocity parallel flow 
ne entrance of the test section. At this point, а 
шискъртав cylinder was introduced into the flow to generate 
аке velocity defect which was then subjected to a constant 
Meee pressure gradient. After leaving the test section, 
the flow was maintained at a constant area through a length 
νι πρ. Chen diffused and directed into the centrifugal 
mower and exhausted vertically. 

= The Inlet 

The inlet section of the tunnel was composed of the 
NEN ШОС ПОО, the plenum chamber, and the nozzle. Incoming 
air was first contracted into a two-foot square plenum sec- 
Ion "his eontraction was accomplished by the inlet bell- 
mouth constructed from eight-inch diameter aluminum pipe 
halves. The corner seams were filled with a sanded epoxy 
шээс to effect a smooth transition of the air into parallel 
flow. The plenum section consisted of four two-foot square 
by one-foot long steel sections bolted together. Cork gas- 


kets were inserted between the sections to prevent the influx 


l2 





Mia? at the Joints and to allow for undisturbed air flow 
eo nese junetions. A four-inch thick section of quarter- 
ШЕСІ honeycomb was installed to straighten the air flow and 
prevent swirling. Next, two copper screens were inserted 
reduce the ambient turbulence level in the flow. Finally, 
the air passed from the plenum chamber to the test section 
through a two-foot long steel nozzle with an exit area of 
forty-nine square inches (seven inches square), and produced 
Bwecontraetion ratio of 11.76. Static pressure ports were 
οτε on the right side of the nozzle at the inlet and 
exit for determination of the exit dynamic pressure. 
3. The Test Section 

The test section consisted of two parallel, vertical 
side walls separated byt textile and adjustable top and 
bottom walls. The flexible walls were constructed from 
one-eighth inch Plexiglas fastened to the nozzle with piano 
EES. As described in Reference 1, these walls could be 
varied to produce a variety of different pressure gradients, 
but for this experiment they were set to create a constant 
adverse pressure а followed by a constant favorable 
pressure gradient. | 

One vertical wall was provided with thirteen slots 
for inserting data sensors. These slots were nineteen 
inches high, three-eighths of an inch wide, and spaced five 
Mehes apart. A filler plug for each slot was carefully 
hand milled, and polished to prevent disturbing the flow 


along the wall when the slots were not in use. A thin 


що 





aluminum strip was wedged over each filler plug to ensure a 
mone tit. When these plugs were removed to insert probes, 
ШОШ Л алп; сегісп seal was used to prevent an influx of air 
through the slot (see Figure 84). 
The other vertical wall of the test section had three 
mecuanpular access holes along the axial centerline for 
ANS tments and measurements. In addition, as one can see 
in Fig. 85, thirty-three static pressure ports were inserted 
along the axial centerline at one-and-one-half inch inter- 
vals. A manometer bank, shown in Fig. 86, connected to these 
static pressure ports provided a visual indication of the 
uniformity of the pressure distribution, as well as indica- 
uns Lc presence of any regions of separation. Ав ап 
additional indication of ro oli e Tilov separation, Several 
ШЕШЕ ОТ vutts were attached to the walls in critical areas. 
Having produced the desired environment for the 
type of air flow that this study was initiated to investi- 
gate, a turbulent wake velocity defect was introduced into 
mimes cunnel by placing a one-inch diameter cylinder in the 


inlet of the test section. In order to ensure a turbulent 
boundary layer at den from ehe cylinder, a thin Strip 
of abrasive, shown in Fig. 87, was attached to the stagnation 
ame of the cylinder. 
Ц. The Blower 
MiGereex bing che test section, the air passed 


through two seven-inch square, four-foot long steel sections 


into a diffusor and then into the blower. A centrifugal 


14 





Power operating at three thousand revolutions per minute 
was utilized to pull air through the tunnel. The air was 
then exhausted vertically through Y sections of steel pipe 
K contained two dampers for the control of the total 


ieee tow through the tunnel. 


B. DATA SENSORS 
meee General Description 
The data for this study were taken with the three 
different types of probes shown in Fig. 88: A pitot static 
pope: a single normal hot wire probe, and an x-wire hot 
(маг probe A detachable traversing mechanism was used to 
Position the probes vertically, and the values were recorded 
ena or on an x-y plotter. 
The Traverse 
"μις apparatus was equipped with two constant speed 
ος to drive the mechanism across the nineteen-inch 
σαι span of the tunnel. The probe was clamped to a 
ог support on the traverse, and the vertical position 
measured by indexing a reference mark on this mechanism to a 
scale fixed to the frame of the apparatus. This scale was 
readable to one one-hundreth of an inch. A rack-and-pinion 
gearing system drove a potentiometer that was also attached 
to the movable support to produce an electrical position 


signal for the abscissa of the x-y plotter. 


19 





e Pitot-Sstatic Sensor 

A Prandtl-type pitot-static pressure probe Model 
# PBA-12-F-11-KL, built by United Sensors and measuring 
Emenundreths of an inch from stem to tip, was used to obtain 
mic ana local static pressures at one-tenth of an inch 
intervals through the wake. These pressures were read on 
two variable angle, water-filled U-tube manometers referenced 
to atmospheric pressure. This reference was chosen to nee 
duce the fluctuations in the water levels. Aces beds dia 
erence 1, taking P-total minus P-static directly creates 
very large fluctuations that make accurate data taking diffi- 
meme Ву Separating the two readings, the fluctuations in 
 ИПапотпесег were significantly reduced, but not 
Баш пасеа. | 

A thirá manometer was attached to the plenum chamber 
Ee he tunnel, and referenced to the last static pressure 
Port on the test section. As one can see in Fig. 89, all of 
these manometers were equipped with portable plastic sights 
ей restricted one's line of sight to a plane perpendicu- 
Шай со the manometer, thus precluding possible parallax and 
assuring consistent readings. 

Ц. The Normal Wire 

This sensor, shown in Fig. 90, was a single-wire 
anemometer probe, Model # 1150 built by Thermo-Systems Inc. 
(TSI), measuring 4 inches from stem to tip. Again the probe 


wes fixed to the traverse for data collection. The lead 


16 





from the sensor was connected to a 1054A-30 anemometer 
nete with a 1051-2 monitor and power supply, also built 
by Thermo-Systems Inc. (TSI) (See Fig. 91 & 92), which 
ο ουδ the voltage changes to readable data. Next these 
Ба меге fed into a TSI 1015C er that was used to 
separate the alternating and D.C. signals (see Appendix B). 
output Of the correlator was connected to an analog 
Gata processing computer which indicated the mean square 
Ec pe output of the correlator. Finally, the output of 
ве 5 apparatus was connected to the ordinate of the X-Y 
Iber. recalling that a potentiometer on the traverse was 
noted со the abscissa. 
шини трие A=Ware Anemometer Probe 

ШЕБЕРІНЕ» ΕΠΟΣ. a Model # 1055 built by TSI, con- 
sisted of two hot wires oriented at 45° to the axis of the 
probe such that the output signals Of the T corresponded 
Со и + vandu- v respectively (See Appendix B). These 
data were processed in the same manner as the normal wire 
pucexcept for the use of several additional algebraic 


ШИШ cbsons of the 1015C correlator. 


ЕН 





III. EXPERIMENTAL PROCEDURE 


GENERAL DESCRIPTION 

Ne” previous study described in Reference 1 had con- 
r (T the test section to produce a constant, adverse- 
favorable pressure gradient. Since this was the environment 
Pmeoutred for this study, no changes were made in the test 
Section except to ensure a tight seal between the vertical 
Acne movable walls by the addition of strips of clay 
Bene che lines of contact. 

ШОУ ООШ collection of data was carried out in three 
Anc Steps. First, the pitot-static system was used to 
весь the total and static pressure data that were intro- 
Mio the computer program, described in Appendix E, 
δα to calculate the flow parameters. Next, a normal- 
ene anemometer was inserted in the test section to gather 
uMsrmacron on both the mean and fluctuating axial velocities. 
These hot-wire values were converted to velocities by equa- 
ting the maximum voltage at each station to the maximum value 
peche velocity calculated from the pitot-static data. 
Finally, an X-array hot-wire anemometer probe was used to 
вес both mean and fluctuating data for velocities in both 
the axlal and the lateral directions, as well as Reynolds 
stresses. а: Ge experiment a constant wind tunnel 
speed was maintained by ensuring that the total pressure of 


the Plenum chamber, minus the static pressure from the test 


ES 





ο ons was no less than 9.90 nor preater than 10.00 


inches of water on a thirty-degree inclined manometer. 


БІ РТТОТ-5ТАТІС DATA 

In this phase of the experiment all data were recorded 
(to an estimated accuracy of one one-hundreth of an inch) 
шгот water-filled U-Tube manometers connected to the pitot- 
static probe. At each of eleven stations downstream from 
the cylinder the static and total pressures were recorded 
at intervals through the wake of one-tenth of an inch. At 
Rs vacion the probe was moved near the upper limit of 
We traverse and situated in the free stream, then the tra- 
verse was marched across the flow, recording the values of 
total and static pressure Screech асре -The initial 
determination of the freestream was accomplished by observing 
that the manometer becomes almost totally insensitive to 
65 in the vertical direction in the undisturbed flow. 
A change in total pressure of less than three one-hundredths 
Dian inch for three 2 ош н о е Was baken со be 
the edge of the wake defect. | | 

As the traverse approached the center of the wake the 
intensity of the turbulence became so large that accurate 
Beadings were difficult to obtain due to fluctuations in the 
Mater level. However, by observing the water levels for a 
long period of time one could determine an average height 
that is estimated to be accurate to three one-hundredths of 


NC The same criterion was used at the bottom of the 


ща 





flow channel, as at the top, to determine when the 


freestream had been reached. 


C. AXIAL FLOW ANEMOMETER DATA 

The probe used in this phase was a single normal hot- 
wire anemometer designed to produce voltage changes corres- 
шие to velocity fluctuations in the flow. When the 
tungssten wire was exposed to the air in the wind tunnel it 
Ксоова by the flow over it, thus requiring additional 
ase со maintain a constant temperature in the wire. It 
ШИЕ агепг, then, that the velocity in the tunnel ís pro- 
portional to the voltage required to maintain a constant 
memocravure in the hot-wire. The probe is placed ina 
Шог Circuit in the TSI module, and its output linearized 


Mor roduce a linear funetion of velocity of the form: 


ВЕС) (түе 


where, 
,. η. - 
u - mean flow velocity 
= песаси Velocity 
апа 
4 
je ou dT = 0 (III-2) 


Meeretore if the DC output on the correlator is selected, 


the output would be 
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1 tou) 


Berner, since this signal was a strong function of time, 
the data were averaged using a low pass filter with a time 
ου L (T) equal to 0.6 secs. The final output of the 


averager was simply 
екп (III-3) 


and this quantity was plotted versus the traverse position. 
The second type of data for this phase was taken with 

En uSelector switch of the correlator on AC, thus producing 

E put corresponding to the fluctuating portion of the 


flow: 
е" = сц! το 


This signal was also a strong function of time but, because 
ο πε nature of the term u', as expressed in (III-2), the 
mal could not be directly averaged. Therefore, in order 
Bo provide a finite measure, u' was first squared and then 


averaged and plotted versus the traverse position, since 


QU AO 


HH 
— 
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Once a plot of e and e' had been obtained, the velocity 
emma be Calculated by solving equations (III-3) and (ITI-4) 


корет 


x 

ο Ξ----- (III-5) 
Ч 
е 


In order to produce accurate values for u and u' in the 
laboratory environment the equipment had to be carefully 
τε Тог each individual flow. Specific instructions 
r LS procedure can be found in the TSI operational 
Manual but, in general, the first step was to adjust the 
control resistance to balance the bridge circuit. This 
adjustment was accomplished by setting the eontrolsresıstor 
шэг аше that would produce a specified voltage change on 
the monitor. At this point the system was ready to run and 
che next step was to Ser Ene linearized operational mode 
Acne monitor, position the probe in the free stream, and 
maximize the frequency response of the wire by setting the 
inductance to the minimum value that would not cause the 


meneame to oscillate. The stability of the system was 
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Monitored during this SD ON 4 dúal-beam oscilloscope to 
mé NO Oscillations were present in the signal. When 

E obe had been adjusted to the flow a digital voltmeter 
шазнисей со check the multiplier-averager for correct opera- 
foe Ihnen the x-y plotter was set to record the signal on 
EN oropriate scale. 

ο. τος the system had been checked for correct operation, 
it was positioned as close to the upper wall as practicable 
and the pen was dropped on the plotter. The sensor was then 
meer) across tne wake to a point near the bottom wall, and 
шшенреп lifted. After changing the paper in the plotter, 
the same procedure was used to make a traverse back to the 
EBENEN "the test section. A total of four traverses were 
made at each station to record both the mean flow and the 
mean square of the fluctuating flow. Each of these quan- 

n ic was recorded at two different scales in order to 


ИСЛ сасе translation into digital data for the IBM-360. 


р. TWO-DIMENSIONAL ANEMOMETER DATA 

The basic operation and accuracy checks of the x-array 
hot wire are very similar to those of the single normal 
wire. However, some additional complications are created 
by having two separate wires on the probe. The terms are 


mow defined, 





As explained in Appendix B, the outputs of the two wires 


(A & B) were 


(D 
| 


= С, (и + у) 


ев = СВ (u = v) 


E op» a given flow, the values of С, апа Сы cannot be 
SS mined easily. In order to find these constants one 


must have 


EN us study the author used a basic physical character- 
Нэг ОТ ріре flow Со set these constants equal, i.e., that 
the lateral velocity must be zero at the centerline. Recall 


шас, асу = о 


Thus the calue of (ει-εμ) was set equal to zero so the 
Bonstant (C) could be defined as in (III-5). Note that this 
Procedure is equivalent to setting V = 0 at the centerline 


Since, 


ο Е = Са РУ) = Clu - v) = 2Cv 
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e, + ев = НУ) + Clg =v) = 2Cu 
The actual position of the centerline was determined from 
memo lcCOu-Static data taken in phase 1. 

Again, the wires were adjusted to the individual flow at 
each station and the data were taken in the same manner as 
@emeribed in the previous section. However, several more 
traverses were required to get all of the data available 
Мосе x-array hot wire. The settings on the correlator, 
Ando he corresponding flow parameters, are listed below. 

All barred quantities were taken with the correlator on 
EN гс sveraged before plotting. The primed quantities 
were taken with hs вовне tana. squared, averaged, 


and then plotted. у 


SETTING RESULT SETTING RESULT 
(A+B) u (A+B) '“ Шы 
(А-В) и, ` 21. v? 

K T pr? (utv)'* 
B πε Bre (u-v)'* 
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TV. DISCUSSION OF RESULTS 


Since this study was intended to carry out an investi- 
EOD of free turbulent mixing in a constant pressure 
gradient, it would be appropriate to discuss how close the 
actual experimental environment approached that desired for 
ши изсийу. It can be easily seen in Fig. 2l that a wake 
MemocitLy defect was initially created by the cylinder as 
desired, and the magnitude of NE in Fig. 41 clearly shows 
шіп спе flow was, in fact, turbulent. However, it can 
also be seen that the adverse pressure gradient encourages 
 Тогпасіоп of a large boundary layer which, as it mixes 
r Stream, interaets with the wake defect, tending to 
Шээс some doubt that free turbulent mixing is still present 
EStpream of Station 6. It is especially significant to 
meme that the turbulence near the wall builds up from a 
negligible value at Station 1 (see Fig. 41-50) to a value 
greater than that induced into the wake by the cylinder 
midway down the test section. As expected; the wake defect 
E ша у begins to mix out'in the favorable pressure 
gradient. 

By comparing the stress calculated from the pitot-static 
data with that from the hot-wire anemometer, one can easily 
ЕБЕ (һас іп both cases the stress is zero on the centerline, 
Men builds up to a maximum as the velocity gradient nears 


the region of its maximum. At the later stations, both 
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types of data show a shear stress sign reversal, and its 
118 begins to increase as the traverse approaches the 
edge of the test section. This trend in the data was prob- 
ШЕ гегіесбіпр the increasing interaction with the boundary 
Meyer as the air flowed downstream. 

Шээс une greater capabilities of the hot-wire probe, 
it was used to take data very near the wall, and the velocity 
De. seen to fall rapidly in the vicinity of the boundary 
layer. The hot-wire data (Figures 1-80) were then used to 
i are Table I to display the value of the correlation 
coefficient 2) and n in the turbulent flow field. Note 
ιν The correlation coefficient oe approaches zero at the 
centerline, and reaches its maximum where the shear is a 
Ens um as predicted by Schlicting (See Fig. 18.4 Ref. 2). 

Considering the relation between kinetic energy and the 
Shear stress, it is easy to see in Table II that both param- 
eters reach a maximum at nearly the same point in the flow. 
uD ring the anemometer with the pitot-static data in 
Appendix F, the maximum shear there απο οσον ас спе 


lateral position corresponding to the kinetic energy maxi- 
ШЕШ Thus, it is clear SENSE least some significant func- 
tional relationship exists between them. The ideal relation- 
Ship would be that postulated by the theory of eddy viscosity, 
that the two parameters differ only by a constant that is 


lotion of x. The slight displacement of the two maxi- 


mums in this case would imply that the relationship between 
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k and C. шэг гт more complex function. Ав ап 
Elec ional attempt to clarify the relation between kinetic 
energy and shear stress, equation (7) was approximated by 
Eme the characteristic velocity as the square root of 


meen rinetic energy, and 


Sz C4 Е) (IV-1) 


AS a result, the relation for Reynolds stress becomes, 


(IV-2) 


The values for equation (IV-2) were calculated and plotted 
EUS traverse position for several stations from the data 
in Table II,which was derived from Figures 21-80. The 
Mestlts of these plots are shown in Figures 93-97. Once 
Ом леге is a definite similarity between the results 
of the approximating equation, (7), and the calculation of 
u'v’ from the hot wire data. In this case the shapes of 
E cacurves were exceptionally similar, but the maximums 
Mere Still displaced by a small amount, preventing direct 
ewn elation by multiplication by a constant. The author 
feels that a more comprehensive analysis of the data in 
this report is warranted by the obvious visual relationship 
of the curves. An in-depth analysis might well uncover 

a Significantly better model for calculating turbulent 


Aw fields than are currently in use. 
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Ana Sheer Stress Parameters 
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The Test Section During Phase 1 
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Station 2 Turbulent Stress Comparison 
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APPENDIX A 
DEVELOPMENT OF THE GOVERNING EQUATIONS 


OF MOTION FOR TURBULENT INCOMPRESSIBLE FLOW 


ШИ спе totiowing equations are the basic Navier- 
Mores equation for incompressible Newtonian flow: 
апоу: 


O (ei) 


momentum: 


гого. у= ү E) 


Ending the DEL terms they are 


2 2 
e + us + узи. + woe rn ас + ο шоо Е 2-3) (А-1) 

: Р ðX ΟΥ 92 
ὃν κ μὸν (ὃν, ὃν 13», 2 Y, 3“ 437v) (дБ) 
9 Ὁ 9х ΟΥ 2, р ду x 3z* 

1 ЭР ο, δ΄ , δ΄ 

Оо 21 11011:5.00151571-0 09) 
3t эх ду dz P эх ду 97 


141 





Now substitute equation (A-1) into (A-4), (A-5), and (4-6) 


Бон сет 
2 2 2 
gu д 2 д 9 TEOT aum дъ 9 w 
w+ = U + = UV + цн = -=> + v(— + — + ) 
ot dx ΟΥ 9Z о 9X x 3y 242 
2 2 5 
= + s uv t 5 v* t 5 vu 22 ас ο ο τις 2-8 ΠΟΤ 
y р ох ΟΥ 92 
2 2 2 
QW 9 9 ο . 1 ӘР 0 W ди д W 
тр + — uw + = ww + Zw” = = ع‎ + v( шаг) 
ot 9х ду 92 р 92 TC 3y* 222 


With p and v constant and using the following definitions 


equation (A-3) becomes 


9 11 ! 9 17 ! 9 u 1 - = 
y CU + u!) + Y + y!) + 210 + и!) 0 (А-7) 


Since the integral over time of any fluctuating quantity is 


zero, taking the time average of equation (A-7) results in 


за ὃν δη. 8 
2156 E ΟΥ ш 0 (A-8) 


дц! + ον + ди! = 0 (A-9) 





substitute the definitions above into equations (A-4), 
then 
с atu + Š (ima po (tu) (avr) + 2 (uru!)(w+w'!) 
ot 9х ΟΥ 92 


ду 
(А-10) 


2 
ООО КЕ un ο ent) σέ 
X x^ 2 922 


1 

| 
TES 
а2| 2 


and after taking the time average 




















ЖЕ) 8-2,9 ,2,9--, 3 9 — , 9 
RI 2 M. t + — БЕЛА t OES + — ty, E: 
St zy + PAL r ү + aye у zu эЧ М (A-11) 
= 2— 2— 2-- 
E 9Р + Орог + 9 u d 3 u] 
p ox эх? a 222 
ema similar manner, one can get 
av ВЕ == πι + 9 Жа E a ye + ο στο DEN zo 
"t + PAL + 3X V S ay" τον ш (А-12) 
me OP 32у, 92у 92у 
-= — + vl + 7722) + = | 
d 9х 9 ο- 
BUE art + vw 2 Š vne a ὃ 2° + Š w: (A-13) 
9% ох 9х ΟΥ ду 92 97 


In the special case for a two-dimensional mean flow with a 


steady mean, equations (A-11) and (A-12) can be written 
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au“ + τν ἂν - -5 © (р+ри' ©) - $, (0 v7) Е ox + = 
(4-11) 

2 212 - 8 (итут) + м3 d ES 

jx Зу p Зу 9х эх? зу? 
(A-15) 


Then applying the standard boundary layer assumptions 


ОЩЕ ща 
Р = Р(х) from aa О 0 BE (21-00) 


and in addition make the following approximations 





p ut“<< P ΠΟ 
N 

ν[5-5 + 2-1 << 2 852 (21-10) 
ох ΟΥ y | 


When these assumptions are substituted into (A-1545), the 


Porlowing equation results: 
a a ёс А а (А-19) 


"Which can be manipulated into the form in equation (A-20) 


meane continuity 


πα τοι, : o EV (A-20) 
dx ΟΥ 9 
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applying the Boussinesq laminarization, the turbulent shear 


becomes 





EE τπτ, cr Eu 
US ; T 5 Cu, (x) 5 y | (A-21) 


EET - -pu'v'!, is the Reynolds Stress, and u. (x), the 
meeye Viscosity, is a function of the flow field. If (A-16) 


(1-21) аге substituted into (A-20) the result is 


5S ML y “Ὁ. ν 24 (A-22) 
2 y еще d зу? 


AN alternate approach to solving these equations is 
possible without the boundary layer assumptions (A-16, 17, 


& 18). Since there are three equations, (A--8), (A-14), and 


—m 0300000000 m e 


ES) and six unknowns u, V, D, ΠΠ. vie. and u'v', three 
ME” ional equations are required. By subtracting the mean 
equations from the full equations as shown in (A-10), three 
additional momentum equations for the fluctuating portion 
πε flow result; however, several additional unknowns 


also appear: 





D ES disp up 
One might also formulate additional equations for these 


variables, but higher order correlation unknowns would 


appear in them, and additional equations would be required 
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ad infinitum. It is customary at this point then to approxi- 


mate vit, ne πο αυ Some function of the mean flow 


solve the system of equations. 
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APPENDIX B 


DEVELOPMENT OF THE EQUATIONS FOR 


THE HOT-WIRE ANEMOMETER 


I. Single Normal Hot-Wire 
Meese assume the following definitions 
e = e + e! 


u Uu + u 


в. since the linearized output of the hot-wire is 


Eooportional to u 


= ш (B-1) 


эг 


οι οσα ου) (B-2) 


BeuLving averages equal to averages, since they must be 


magependent of the fluctuations, 


(B=3) 


Ф| 
il 
c 
= 


Æ Solving for u 


(8-4) 


Z 
1 
оғ 
Ф| 
| 
το 
(D 
Ф| 


Ф 
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ana 


NSR 
I 
Ф|19| 
| 
T 
wn 


Now setting the fluctuations equal 


a y? (8-0) 
u! z 5 е! 
NOD 

ος 


which from (B-4) becomes 


2 
ЕН τ. 
ше. -6% οὐ" (B-7) 
Е 
oT 
ut^ etc 
ше) 5 1-2 en 
Ve Ce 


II.  X-array Hot-Wire 


Now assume in addition, that 


B 
(B-9) 
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then since the wires are separated by 90° 








ο Ισ EOS 0 


tu 
— > — T. 
4 
L ucos® + vsin06 
ο + V 
| 
wire A wire B 


Bemznote that Ө = 45° so 


е, = С, ( u + v) 


A 
(B-10) 
eB = Cp (u-v) 
At the centerline of the flow the lateral velocity is 
Epeoximately zero, so ед + en and, therefore, σι = Св = С. 
Нэг this reason, once the probe is adjusted to force ед 5 ер 
at the centerline the following relations hold 
ед + C (u * v) 
(B-11) 
ев = (00220) 
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But if these equations are summed and solved for u 


E d. 
u = 56 le, + ep] (ВЕЩ 
emo Similarly 
ee А Go) 
20 p B 


Then if B = 1/2 C, and the definitions of the parameters are 
substituted into (B-12) & (B-13), the resulting equations 


are 


сет = доо ' 2012 ' T 
u + u Ble, + е, + en + en ] (В-14) 
ү "= En Ни үл 2 ! = 
Ev Ben te, ep = €p' J (BEIS) 
шити В-Т, set averages equal to averages and 
ШІ = Ble, + ев] = Ble, + ep] (B-16) 
NAS Ble, ын ев] = Ble, ES en} (B-17) 
Bine fluctuations equal to fluctuations, 
== — ! = 
m ο ο ο IS ELE, + е)" | (В-18) 
= 0-2 ' = = i -- 
v' = Ble,? - eg'] = Blle, - eg) (B-19) 
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Then square and 


с 
И 


< 
! 


Returning 
primes 
en = C(u' 
ен! = C(u! 
Amare and average 


+ y! 


average to get 
дар, па 
8 [(e, t ey)' 1] 


8^[(e, - ey)'^] 


) 


- y!) 








then subtract 





"o 


SA 


or 





(u'^ 4 2ü'v* + у! 


=) 


wen mama =e mm с ч 


(ц'“ A лд 








5 (ед 


Now add equations (B-23) 


a] 


— 


(B-20) 


(B-21) 


to equations (B-11) and setting primes equal to 


(B-22) 


(en 


(B-24) 


(В-25) 





and 





or 


equating (B-23) 


(ед 
Ипеге 
"Te - 
or 
> [To 
so 
Ве = 














en!” = әсери: с + E 
12 т 2 ща ща 
V 28 IE} + en 1 (8-23) 
and (B-21) 
по ВЕЕР ПОР = 2 
V B (е, + ед)! 4 den ең)! 1 (В-24) 
and (B-24), and rearranging 
2[e,'^ + е! 41 
Sl (B-25) 
2 2 
1 x 1 
+ ев) + (ει ев) 
ze 2 > 
ВА = 1/2plu, + VE ] 
ша 1.2 1.2 5 
Е ер en). 
р 


- Р 


2 ("те е) 


p[ (e, + ο, + (е, - SA 





and 


2 (те = a 


is a known array of data from the pitot-static phase of the 


EN y . 
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ANDA С 
DEVELOPMENT OF AN INTEGRATED FORM OF ΤΗΕ, 
MOMENTUM EQUATION FOR THE CALCULATION OF 


THE LOCAL TURBULENT SHEAR COEFFICIENT 


peso. the basic equations of motion for a two-dimensional 


mean flow are 


continuity: 
u, + Vy = 0 (C-1) 
momentum: 
1 1 
toy == АРЕ + — 0-2 
uu. u, e E ( ) 


X y 
y y 
f u, dy = =f Vv ау 
0 Олар 
y y 
fu. dy = - [v + f(x)] 
x 
0 
thus 
y 
v(y) = - Р u. dy + poe 
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But v IS Zero. SO 
1 i 


y 


Же е u, ду (С-3) 
0 


Now simply and multiply equation (C-1) by u and add to 


(C-2) to get 


T 
1. E „у (22 
Х В y p < р 


where the terms on the left can be recognized as differentials 


of the form 


u 2 
anu = (u )x 


τ у = 
vu, u а (чу) 


thus equation (C-4) becomes 


i “у 
2uu, + (uv) = - ЕР + —- 
X y P X P 
since V(0) = 0, when this equation is integrated from 
w= 0 to an arbitrary point y = ve 
y a 2 
Iu uc v7 dy 2 2uu. dy + uly_)v(y_) 
0 X 0 у 0 χ т D 


and the result, when equation (C-3) is substituted for 


Ey, is the following relation 
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Ут Ут US s 
КИ ο ο (у, ау = = Г —dy + 3 
0 X τρ X 0 p 





ЕП сап be rearranged into the following form, 


y y 


i 1 Т 1 
! u,L2u-u (y, )]ady pus нг. 5 Dos) (0-5) 


0 


And can be solved for the coefficient of shear at У 





2 
(С. = 5 т) 
PU a ` 


y y 
HER E 2 Б 
C y.) = 5 2 ВЕ СИ uy L2u-u(y. ) Jày] 


pu u 0 
e e (0-6) 








An alternate approach to this development would be to 


substitute equation (C-3) into (C-2) 


e y 
Qu Qu Qu 


ат 
ον ο) (ὃν (0-1) 


сін 
Q? i Q2 
AET 
+ 
о |ы 
сс 


Then non-dimensionalize by division by the freestream 
Néelocity squared, and use the product rule of the calculus 


to get the following form of the momentum equation: 





нв. + (Uy Чех ща А) 17 aya) (EX) A 
NIS. ^  ΄ зто ет о 
1 C. (C-8) 
= 5 Р + ف‎ 
χ 
pru 2 


с», 





K re с = 21 m сг 10502105 Х, and y indicate 
T 2 


D 





Em c! differentistion with respect to x, and y respectively. 
Now, if the condition that equation (C-8) be satisfied at 
роит іп the flow is relaxed to allow it to be satis- 
feed оп the average between y = 0 and an arbitrary point 

Y then (C-8) can be integrated by parts to get the 


w Owing relation Гог = 





EECy ) Ут uly_) u 
Т u T u u ex 
————— = Ј [2ίτπε)------]ζ(--) +(—) (——) Jdy 
2 0 Че На БЕ На Че 
(0-9) 
y 
Тт 
+ f 1 5 Р, αν 
0 рч 


В: the original definition of the shear term 





Memsubstituted back into equations (C-5) and (C-9), the 


following relations result: 











y y 
Шу е f u, Lu (y ) - 2uJdy - = P. dy (C-10) 
T p X 
0 0 
y 
шиг. u 
ШҮ! = и ^ / Т--2(2-10(2-) 4622-1629 Jay 
0 ας e ex e e 
ee (C-11) 
DE ENS dy 
0 ρ 


1: 





Party, го check the validity of these equations, let 


p. Evo infinity in the limit and assume that 


Then, equation (C-5) becomes 


200 uu, dy τα. / u. dy - f ας U, ду = 0 
0 D 0 X 


2 


meee vie shear, and therefore 6 CS το ο io tana” μας 


00 2 | со со 
эг ОК М idy = гоа ау=о 
0 : ps 0 Нэг: 5 0 = 
эг 
г Гад) - (u u) + uu τ. αν = O 
0 x ας e, e x 
us, combining terms 
2 u u и u 
a) Cy i) UU dy - 0 
0 πμ -- X 0 E NEG 
(C-12) 
and if the flowing parameters are defined 
mu и x et 
6= f —(—-1)dy; 6 = f (— - 1) dy 
e e 0 e 
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then equation (C-12) may be written 


[u, 9] + uu ὃς - 0 
χ X 


whieh when the first term is differentiated becomes 


27008 А ы 
о 
x > 
E... 2 

ша ту 13 пе Бу и 
ae 

+= [20 + 6%] =0 (C-13) 
X ц 


which is simply the well-known momentum integral equation 
NS boundary layer. A similar procedure will reduce 


eauation (C-9) to (C-13) as well. 
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APPENDIX D 


DATA CHECKS 


moe experimental study, it is very important to have 
Dime idea of the quality of the data that results from the 
ΜΝ. rilment. In this Appendix, several methods of checking 


the accuracy of the data in this study are discussed. 


CONTINUITY CHECKS 


me the continuity equation, it is apparent that, 
Since the mass flow must be constant, 
wall 

í udy = constant (D-1) 

0 
ээ гс ї6 equation was added to the program described 
uu Dpendix E. The results of this check for the pitot- 
static data varied as much as 30% over the entire test 
0, yet, if the largest and smallest value are neglected 
the range is narrowed to + 10%. The author believes that 
meen Of this difference is due to the inability to position 
bEnecpitot-static probe near the wall in some cases, while 
Mothers, the head of the probe could be placed adjacent 
(пе wall surface. This problem is especially critical 
at Station 1 where the boundary layer is thin, and mass 


и епеар (пе wall, at free stream velocity, comprises а 
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significant part of the total mass flow in the tunnel. 


"Πο computer output in Appendix F. 
MOMENTUM CHECKS 


Recall the integral momentum equation derived in 
Appendix С, 


й u 
e 
go + — [20 + 6J = 0 (06-13) 
u 
e 


which is a special case for the. free stream where C is 


meow taking this equation as 


глад 
El = ДЕ 
du 
ERU e 
Ба = а dx 
e 


ντ ου + ὦ 
then the non-dimensional residual of this equation was 


Programmed as 
EI - (ЕШ Е E2 % ЕЗ)/(Е1 - Е2 % ЕЗ) (Da2) 


and its values are listed in the computer output. 
Tri general, this ratio showed the data to correlate 
fairly Пеши ύπε restrictions imposed by the momentum 


integral equation with an average residual of less than 
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25% and as little as +15% for most stations; however, 
τος 7 data was completely out of line with a residual 
greater than two. The primary reason for the failure of 
Eon 7 to check satisfactorily was an inclusion of a 
meat icant portion of data at the upper portion of the 
traverse that reflected boundary layer interference. The 
Sion of the resulting velocity decay in the calcula- 
шой ог the momentum, and displacement thicknesses created 
шини се еггог in their values. This error is almost 
certainly ПОШ И 5с ους 11 Е:0101У high values of Station 8 


as well. 
ENERGY CHECKS 


Ea calculated the local turbulent shear coefficient, 
Mas possible to also integrate to the residual in the 


energy integral equation, 


=2 ЛЕ (¿May 22 


Dre the energy thickness 93 is defined as 


со 


6, = 2 u(u, ^ z u*)dy (D-4) 


Therefore taking the term on the left as A, and that on 


the right as B, the relation was programmed and, 
ZETA = (A - B)/(A * B) (D-4) 
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mer results of this data check were inconclusive. The 
magnitude of one of two terms above (A and B) was so much 
Mareer tnan the other that the difference over the sum was 
οιοντο an error in the programming vice an inaccuracy 


in the raw data. 
HOT-WIRE CHECKS 


In Appendix B it was shown in equation (B-25) that these 
Bons were set equal to n, and values for n are listed in 
Нэг 1. It was immediately obvious in this case that the 
mo although not large, was directly related to the 
Meet of the turbulence in the flow. Thus, the magni- 
MO: the lack of algebraic correlation in this check 
may have been due to an inability of the equipment to 
handle the rapid variations caused by severe turbulence. 
IS important to note here, that a change of as little 
as 10% in the individual values of the parameters used 
to calculate n would bring even the worst value of n back 
to within 10% of 1.0. The author felt that the data from 
Phase III of the experiment was well within the bounds of 


acceptable experimental data spread. 
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NELENDIX E 
DATA REDUCTION COMPUTER PROGRAM 


A. GENERAL 

ши Fortran TV computer program included in this Appen- 
EN s a revised form of that developed by the study of 
Reference 1. The program inputs the pitot-static data of 
ШОО О of the experiment, and a limited amount of data 
from phase II to calculate the wake momentum thickness, 

ФИ асепепс thickness, energy thickness, local turbulent 
E coefficient, and Eddy Viscosity. 

' The inception of the present study the portion of 
the program dealing ᾿,''ν, πε ест соев тастен 
που produce accurate values for С... ines author of 
Memerence 1 was unable to find a reason for this error 
mee the equations appeared correct, and there were по 
Obvious coding errors. He postulated that the unreasonable 
Values were a result of an inability Cf aenorweard finice 
Gifference scheme to accuravely approximate the partial 
derivatives in equation (0-9) Gf ПГ μα. ο sce this 
appeared to be a reasonable theory, this author attempted 
mo, Obtain values for the derivatives by two additional 
methods: 

1. A least squares polynomial approximation was fitted 


to the data at a given position y from the 
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ΤΠΕ cuen this polynomial was differentiated. 

2. A central difference scheme. 

Unfortunately, these improvements resulted in no signi- 
Mano” change in the values of t. Thus, the values for 
E various derivatives were printed out in detail, and 
ЕП order of magnitude analysis was carried out. At this 
Meant, it was not difficult to see that the values of the 
 шесет Ч were exceptionally smooth in the x-direction, 
and obviously not a problem to approximate, yet the values 
calculated by the IBM-360 were not even similar to those 
calculated by hand. The author then attempted to follow 
Mier instruetions programmed for the computer and quickly 
ED ovepred that a subtle error was hidden in the program. 
Ecce each station had a different number of points, the 
Ex corresponding to a point y above the centerline was 
EU rerent at each station; therefore the arrays requiring 
differentiation were read into a dummy variable where 
NEED corresponded to the centerline for each station. 
ШИЕ correction produced values for the local shear coeffi- 
ИШЕТ that were of ІШІНЕ right order τος, and the 
Nentral difference scheme was considered to produce the 
best values for the derivative. Some of the hot-wire data 
Was also analyzed using this program, but the values for 
C. were not very similar to those calculated from the pitot- 
static data. Since the pressure term was found to be the 


driving force of the equation for Cx during the order of 
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magnitude analysis, the author felt that approximating 


ШИЕ static pressure data with the data from phase one would 


EEunlikely to produce realistic values for the parameters, 


and no additional attempt was made to use the IBM-360 to 


reduce the hot-wire data. 


B. NOMENCLATURE 


The following parameters and arrays were used in the 


iam; their respective dimensions are shown in 


parenthesis: 

ALV ο ο used to store interpolated values of 
dynamic pressures (feet) 

AMIXL Jocal Prandtl mixing length in the wake (fee) 

CTAU ioc turbulent Shear coefficient for each point 
in the wake 

DELRAT maveo Of SUCceeding DELSTR's to that at the first 
data station 

DELSTR wake displacement thickness (feet) 

DELX Seema α.στοπςος between succeeding probe tip 
positions (feet) 

DIAMTR diameter of wake generator (inches) 

EDDY eddy viscosity 

EPSIL residual in momentum integral equation when 
experimental values for THETA and SHAPE are 
Subst реа into it 

LABEL label for subroutine DRAW 

MN number of data points laterally through the wake 

MX number of axial data positions 

MY number of points (always odd) into which the 


wake is divided 
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PATM 


Q 
SHAPE 
ШНЕТА 


HET AP 


THETAV 


ШНКАТ 


BILE 
TITLE! 


U 


URATIO 


XIN 


XOUT 


XPOS 


FOUT 


LPS 


YQ 


MU(X) 


MM 


array of local static pressures for each point 
in the wake (psf) 


atmospheric pressure (inches of mercury) 


freestream dynamic pressure (psf) 


wake shape factor (DELSTR/THETA) 


wake momentum thickness (feet) 


inclination of static pressure manometer 
(degrees) 


inclination of dynamic pressure manometer 
(degrees) 


ratio 


SO ce cane TER IAS to that at the first 


data station 


title for subroutine DRAW 

τοις for subroutine DRAW 

array of velocities for each point in the wake 
(feet per second) 

5 οἱ local wake velocity to local freestream 
velocity 

array of probe traverse positions through the 


wake (inches) 


interpolated probe traverse positions for each 
of the MY points in the wake (inches) 


probe 


tip position measured from downstream 


side of the wake generator (inches) 


interpolated manometer readings (for YQ and YPS) 
Шо есеп о: the MY points through the wake 
(inches ) 


Static pressure manometer slant distance height 
difference (inches) 


dynamic pressure manometer slant distance 
height difference (inches) 


kinematic viscosity 


index of the maximum velocity 


πο 





MT 


j 


TRATIO 


ADP 
D 

NUM 
MBM 


MAD 


MID 


index of the maximum velocity in the upper 
του οὗ πας 


dummy 
dummy 
dummy 
dummy 
array 


array 


array 
array 
array 
array 
of MN 


of MM 


data 

Cie Or differential 

of URATIO O 1 1 
oi Tor differential 


oT поло О ШОГО Eddy Viscosity 


array of MID 
array of MT 


index of centerline 


С. SUBROUTINES 


INTEG 


QTFE 


QSF 


FILL 


PREIN 


ШОР 


TIT 


D. USAGE 


uses Simpson's rule to numerically integrate the 
Ши гс го Ве it must be supplied with K points 
(always odd) equally spaced, DELX units apart 


рога Вороне ror integration using 
trapazoidal rule 


Doo library subroutine Tor integration using 
Simpson's Rule 


used to fill in arrays to maximum dimension with 
freestream values 


determines midpoint of asignally dimensioned 
array | 


determines MT, & MM 
takes the entire array and puts the values above 


the centerline into a second array in the 
Шил parameters 


The details on the use of this program are listed in the 


comments at the head of the program. 


The deck requires 120K 


to compile and 230K to run. 
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